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r^ ■ I. INTRODUCTION 

(N , 

Up to now most of the hadrons observed experimetally can be interpreted as the qq/qqq states in the quark 

T^ , model [l|, [2|. However there has accumulated some evidence of the exotic state with J^'~^ = 1 ^ [SJ-Q- Such a 

^p^' quantum number is not accessible for a conventional meson composed of a pair of quark and anti-quark in the non- 

O ' relativistic quark model. Sometimes these states are named as exotic states although all the J^^ quantum numbers 

QJ • are allowed in QCD. 

T^ ' For a neutral quark model qq state, J = ensures L = S hence C = (— )^^'^ = +1. Therefore, there exist two 

possible exotic states with J^*^ ~ Q~~ and 0^ . It's also interesting to note that the J^*^ quantum number of the 

^\i ' local operators composed of a pair of the gluon field strength tensor is cither O'*'^ or '". 

^ , On the other hand, the tetraquark operators may carry the and 0^* quantum numbers. In fact, the state 

On ■ was investigated systematically using the tetraquark currents with the QCD sum rule method |6|, |7[. As a byproduct, 
CTN , it was noted that there does not exist any tetraquark interpolating current without derivative for the J^*^ = 0"' 
case. 

With the similar formalism, one may construct the possible 0"' tetraquark current by introducing derivatives. 
C^ . There are two kinds of constructions either with the qq basis or with the qq basis: {qq){qq) and {qq){qq)- However, 
^— ^ ' they can be related to each other by the Fierz transformation [y|. In this work we use the first set. With these 
^ , , independent 0^ currents, we perform the QCD sum rule analysis and extract the masses of the corresponding 
currents. 

This paper is organized as follows. In Sec. [Hi we construct the tetraquark currents with J^'-^ = 0^ using the 
diquark {qq) and antidiquark (qq) fields. In Sec. IIIIl we calculate the correlation functions and spectral densities of 
the interpolating currents and collect them in the Appendix |B] Wc perform the numerical analysis and extract the 
C^ ■ masses in Sec. II VI for the light and heavy systems respectively. The last section is a brief summary. 

II. TETRAQUARK INTERPOLATING CURRENTS 

It was shown that the J^"^ = 0"' tetraquark interpolating currents without derivatives do not exist (6|. So in this 
work we construct the 0"' currents with the derivatives following the similar steps as in Ref. Q. We first construct 
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two independent tetraquark fields: 

Aabcd = {qLC7'"q2b){q3c D^ CqJa) , (1) 

Pabcd = (gfa C7'' 75926) (gSc D^ IbCqJd)- 

where qi_4 represents the flavor of quarks, and a — d stands for the color indiees, Dfj.~ Dp, — -t'p, Dp, = 5 ^ + igA°jt°-. 

It is understood that the index c is the color index of I qD ^ j . In Eqs.([T]) and ^ we have used the shorthand notation 

to simply the expression. 

To compose the color singlet tetraquark currents, the diquark and antidiquark should have the same color and spin 
symmetries. Therefore the color structure of the tetraquark is either 6 ® 6 or 3 (8) 3, which is denoted by labels 6 and 
3 respectively. Details can be found in Ref. Q. Considering both the color and Lorentz structures, we can obtain the 
currents with J^'^ = 0^ : 



?7i(x) = qiaC-f^''q2b{qia Dp Cq2b + qib Dp Cq2a) - qia^ Dp q2b{qial^''Cq2b + qibl^Cq2a), 
V2{x) = qlaCY''q2b{qia Dp CqJb - qib Dp Cq2a) - qlaC Dp q2b{qiaY''Cqlb ~ qul'^CqJa), 

m{x) = qla'^i^i5q2b{qia Dp i^Cqlb + qu Dp ibCql^) - gf^C Dp i5q2b{qiai^i5Cqlb + qui^ibCql^), 
T]4{x) = qlaCY''ibq2b{qia Dp ibCq^b ^ qib Dp ibCqla) - qicfi Dp i^q2b{qiai^ibCqlb - qibi^izCqla)- 



(2) 



III. QCD SUM RULE 

Consider the two-point correlation function in the framework of QCD sum rule 

n(g2) = J /d4a;e*?^(0|r7/(a;)77^(0)|0), (3) 



where 77 is an interpolating current. At the hadron level, the correlation function n(|3^) is expressed via the dispersion 
relation: 

^{p') = r ^y . ds, (4) 

J a s-p- - le 
where 

p{s)^Y.^{s-mlmTi\n){nW\0) 

(5) 
=fx5{s — m-^) + continuum , 

where mx is the mass of the resonance X and fx is the decay constant of the meson: 

(OhiX) = fx. (6) 

The correlation function can also be calculated at the quark-gluon level using the QCD operator product expansion 
(OPE) method. It is convenient to evaluate the Wilson coefficient in the coordinate space for the light quark systems 
and in the momentum space for the heavy quark systems respectively. In our calculation we consider the first order 
perturbative and various condensates contributions. In order to calculate the gluonic condensate, it is convenient to 
work in the fixed-point gauge. The massive quark propagator iS{x, y) in an external field in the fixed-point gauge 
is listed in Appendix A. The quark lines attached with gluon contain terms proportional to y, which we can ignore 
in the current without derivatives. We keep these terms throughout the evaluation and let y go to zero only after 
finishing the derivatives. The n(p^) can be written as: 

foo OPE( \ 

U^'-'^ip') = / ds P ^ ^'\ , (7) 

where the TOi and 7712 are the mass of the quark gi and (72 respectively. In order to suppress the higher state 
contributions, we perform the Borel transformation to the correlation function, which improves the convergence of 
the OPE scries. With the quark- hadron duality, we obtain: 

/•oo 

n(M|) = /ie-™-/^^- = / dse-^/^^-pO^^(5), (8) 



where Sq is the threshold parameter and Mb is the Borel parameter. We can extract the meson mass mx'- 



a(Tfe)n(M|) _ /4T,„,+„,).rf-"^/^^S^p(.s) 



'x 



n(M|; 






(9) 



For all the tetraquark currents m Eq.([2]), we collect the spectral density p'^^^{s) m the Appendix. The quark 
condensate {qq) vanishes due to the special Lorenz structures of the currents. For the q ~ u,d, we do the calculation 
in the chiral limit ruq = 0. Since the contribution of the three gluon condensate {g^fGGG) is very small, we consider 
only the power corrections from the following condensates {g^GG), {qgcr-Gq), {qq^ and {qg(T-Gq){qq). We list several 
typical Fcynman diagrams in the Fig. [T] 
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FIG. 1: Some typical Feynman diagrams of the correlation fuctions. 



IV. NUMERICAL ANALYSIS 



In the QCD sum rule analysis,we use the following values of the quark masses, coupling constant and various 
condensates dH^ISl: 

mc{mc) = (1.23 ±0.09) GeV, 
rribimb) = (4.20 ± 0.07) GeV, 
{qq) ==-(0.23 ±0.03)3 Q^yS^ 
{qga-Gq)^-M^{qq), 
Ml = (0.8 ± 0.2) GeV, 
(ss)/(g(z)=0.8±0.2, 
{g^GG) = (0.88 ±0.13) GeV^, 
a^l.TGeV) = 0.328 ± 0.03 ± 0.025. 

The Borel mass Mb and the threshold value sq are two pivotal parameters. The working region of the Borel mass 
is determined by the convergence of the OPE and the pole contribution. The requirement of the convergence of the 
OPE determines the lower bound MBmin of the Borel mass, and the pole contribution determines the upper bound 

In this work, there is no contribution from the quark condensate (qq). The correction from the condensate (qgcr ■ 
Gq){qq) is the most important numerically. Its contribution is bigger than that from the gluon condensate (g'^GG), 
mixed condensate (qga ■ Gq) and the four quark condensate (qq)'^. The mixed condensate is also very important 
numerically for the currents 771.3. We require that the condensate {qgcr ■ Gq){qq) be less than one ninth of the 
perturbative term to ensure the convergence of OPE, which leads to the lower bound of the Borel parameter working 
window. 
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TABLE I: The threshold values, Borel window, Borel parameter for the different tetraquark currents. 



The pole contribution (PC) is defined as 



PC = 



-'4(7ni+??i2) 



,dse-''/'^'Bp(s) 



L 



4(mi+m2)^ 



dse '^/^^Bp(^s) 



(11) 



which depends on both the Borel mass Mb and the threshold value sq. sq is chosen around the region where the 
variation of mx with Mb is minimum. Requiring the PC be larger than 30% '^ 50%, we get the upper bound MBmax 
of the Borel mass A'Ib- We list the working region of the Borel parameter for the four currents with different quark 
composition in Table HI For the r}'l_^, we get the upper bound of Borel parameter Mb by requiring the PC be larger 
than 30%. For 77j_4 we require the PC be larger than 40%. The masses are extracted using the threshold values sq 
and Borel parameters M^ listed in Table ID The last column is the pole contribution with the corresponding sq and 
Mb. 

For the light tetraquark systems, there does not exist a working region for the sum rules. Even in the extreme case 
that the pole contribution is ~ 30% and the contribution of the condensate {qga ■ Gq) (qq) is around the leading order 
contribution, the lower bound MBmin is still much larger than the upper bound MBmax- In other words, there is no 
working region for light quark systems. As shown in Fig. (|2l3p . the extracted mass grows monotonically with sq which 
implies the 0^ state does not exist below 2 GeV. We note that the light J^'^ = 0^^ state does not exist either [y]. 
The 0^ and ^ channels are in strong contrast with the 0++ case, where there exist stable tetraquark QCD sum 
rules and the extracted scalar meson masses agree with the experimental scalar spectrum nicely |ll| . 

For the heavy systems, the variation of mx with sq and Mb is presented in Figs. (|4))-([TT|). All the sum rules are 
very stable with reasonable variations of sg a-nd Mb- The presence of the two heavy quarks reduces the kinetic energy 
of the tetraquark system, hence helps to stabilize the sum rules. Numerically, the masses of the 0"' states are slightly 
larger than those of the 0~~ states Q- 



SUMMARY 



The exotic state with J^'-' = 0"' cannot be composed of a pair of quark and anti-quark. In order to explore 
these exotic states, we have constructed four tetraquark interpolating operators. Then we make the operator product 
expansion and extract the spectral density. Because of the special Lorcntz structures of the currents, the quark 
condensate (qq) vanishes. 

For the light tetraquark systems, there does not exist a working region of the Borel parameter and threshold for 
all the derived sum rules. It seems that none of these independent interpolating currents supports a resonant signal 
below 2 GeV, which is consistent with the current experimental data [l|. In contrast, there exist very stable QCD sum 
rules constructed from the tetraquark interpolating operators in the scalar channel. The extracted scalar spectrum 
agrees with the experimental data nicely |ll| . 

For the heavy quark systems, the 0"' tetraquark sum rules are quite stable. The presence of two heavy quarks may 
render the kinetic energy of the tetraquark system, which is helpful in the formation of bound states. The extracted 
masses from the four interpolating currents ?7j„4 are around 4.76 ~ 4.96 GeV for the charmonium-like states. For the 
bottomnium-like 0"' states, their masses are about 11.2 ^ 11.3 GeV. It's very interesting to note that the mass of 
the 0"' charmonium-like state extracted from the tetraquark sum rules is numerically quite close to the mass of the 
0^ hybrid charmonium extracted on the lattice jla[l3j. 

Because of the special "exotic" quantum number, the 0"' charmonium-like state does not decay into a pair of 
particle (H) and anti-particle {H). There are two types of 0"*"" states with different isospin and G-parity: l'~^ = 0~ 
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FIG. 2: The variation of mx with Mb (left) and so (right) for the current r\[. 



1 



- Mj=l.lGeV 

- Mj=1.2GeV 

- Mj=1.3GeV 

- Mj=1.4GeV 



5 7 9 

i„ [GeV^] 



13 



2.0 


v::::v::v::::v::v::::::-:EE^:^ 


% 1.5 

O 








1 '-^ 
0.5 




.t„=5GeV= 

.i„=6GeV- 

i„=7GeV= 




n n 









1.10 1.15 1.20 1.25 1.30 

Mb IGeVJ 



1.35 



1.40 



FIG. 3: The variation of mx with Mb (left) and sq (right) for the current r\\. 
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S-wave 


P-wave 


0" 


Xci(lP)fti(1170)... 


D"(1865)Di(2420)" + c.c.,7:'*(2007)"I>S(2400)" + c.c, 

77,(lS)/ii(1170), J/V(15')/o(600), 

J/^/o(980), J/V^/i(1285), Xco(lP)^(782), Xci(lP)a;(782), 

^(2S)/o(600), V(3770)/o(600) . . . 


1+ 


J/^(lS)7ri(1400), 

J/V'(lS)7ri(1600), 

Xci(lP)fei(1235)... 


D"(1865)I?i (2420)" + c.c, D*(2007)"DS(2400)" -fee, 

D*(2007)°Di(2420)''+ c.c, 

r;c(lS')bi(1235), J/V'(lS')ao(980), J/V'(lS')ai(1260), 

Xco(lP)p(770),Xci(lP)p(770)... 



TABLE II: The possible decay modes of the 0"'" charmonium-like state. 

and I'^ = 1+. Only a few S-wave decay modes are allowed. Some possible two-body decay modes are listed in 
Table HT] Replacing the D meson by B meson, one gets the decay patterns of the bottomonium-like states so long 
as the kinematics allows. The 0"' state may be searched for experimentally at facilities such as Super-B factories, 
PANDE, LHC and RHIC in the future, especially at RHIC and LHC where plenty of charm, anti-charm and light 
quarks are produced. 
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Appendix A: The Momentum Space Propagator 

The fixed-point gauge is defined as: 

ix-xorAlix)=0 (Al) 

wliere xq is an arbitrary point in space wliicli can be cliosen to be the origin. Then the potential A" can be expressed 
in terms of the field strength tensor G^y{G^i, ~ ^G" ) |14l.ll5|: 

A^{x) = / tdtG^^,{tx)x'' 

Jo (A2) 

l-x^G^^iO) + Ix^'x-'Da^G.^iO) + Ix^x^x'^D^DpG^^iS)) + . . . , 
Z 6 o 

Denote the massive quark propagator between tlie position x and y in tlie coordinate space as iS{x,y). The massive 
quark propagator in the momentum space is |14{ : 

iS{p) = iSaip) + iSgip) + iSgg{p) + . . . , (A3) 

where iSo{p) is the free quark propagator: 

tSoip) = -^^, (A4) 

p — m 

where p = 'y'^p^,, iSgijj) is the quark propagator with one gluon leg attached: 

iA" ^„ at^-{p + m) + {p + m)a^- , » A" r2y^Y{p + rn) yf^Y ] .,,^ 

^^aiP) - l^9sG^, (p2_,„2)2 + aT^^^^"^ V {p^-m^f p^^m^\ ^^^' 

iSgg{p) is the quark propagator with two gluon legs attached: 



lbgg[p) ~ -^YY^" ^'P "" (p2 - to2)5 ^■' +*' ^ ■' ' 



l^"^„2^a^6 P + ™ 



'9sGnpG' 



4 2 2"" >'P '"'{p-^-m'^Y 
where ft^'^---°'P = -f'^{p + ra)Y{p + m) . . .j°'{p + m}j^{p + m) 



y^ir" + r") + yT" - ly'y^p' - m^) 



(A6) 



Appendix B: The Spectral Densities 

111 this appendix, we list the spectral densities of the tetraquark interpolating currents. For the light quark systems 
(91792 = u,d)^ the spectral densities are: 

^s^ 17 a {9'GG)s^ {qqfs^ 799 {qg s<J ■ Gq){qq)s 



512007r6^ 108 TT^ 184327r6 6tt^ TGStt^ 

^ s^ T_a {g'GG)s^ _ {qq)\s^ _ 245{qgsa ■ Gq){qq)s 

^^^^' 1024007r6^ 54 TT'^ 184327r6 127r2 7687r2 

^ s' r^a _ {g'GG)s^ _ {gqfs^ _ 7{qgsa ■ Gq){qq)s 

^^^^' 512007r6^ 108 tt'^ 184327r6 Gtt^ 127r2 

^ s" T_a {g^GG)s^ _ {gqfs^ _ Uqg^a ■ Gq){qq)s 



(Bl) 



1024007r6^ 54 tt' 184327r6 127r2 247r2 

For the heavy systems {qi = u,d, 52 = c,d), the spectral densities are: 

p{s) = pP'^^^s) + p<^^>(s) + p'^^^^Hs) + p<99)'(s) + p(9Gg)(g9)(g) (g2) 

For the condensate {qq){qGq), it contains two parts: one part could be written as p, and the other part couldn't, 
which we perform the Borel transformation directly. Therefore 

/•CJO 

For the interpolating current 771: 
384 



"sXij- r 

pP-*(s) = _ [(I6p[i5(s) + mVfi4(s) - 2m2pfi4(s) + 6m'pi,^{s) - pf,,{s) + 2p^,^is) - Ap'^s)) 



a ,17 T 25 2 / / \\ 

+ -[-^Pii5(-V - ^57™ Piu(-V) 



(B4) 



pf''''\s) = i^hpUs) + ^pUs) ~ \p^..^{s) - ^^p^,,,{s) + '-pUs) + \pUs) '-^pU^) 



n'' 



39 5 

- yPL3(s) + 1152mV(34(s) - 20m2p^24(s) - 16m^pi^^{s) - TGSmVLls) " ^Pi23is) 

95 5 

+ 192m2(pf44(s) + pfi4(s)) - 96m^{p?^^is) + p^.^is)) + 1152m* p{^^is) - ^pfial^) + ^m'p^^si^) 

in 5 5 91 

- 192m%p^,,{s) + pfi4(.)) + -m'p^,,{s) - -P2^i4(^) " ^^^'P224{s) " 240mV3^i4(s) + ^PnsC^) 

777^ 2Q SI 8 

- —^^23(5) - yP[24(-^) + 160m2pf34(s) + 3072mV[45(s) + lOm^p^^Js) - -pf24(s) " 3P?i4(5) 

- 192m2pf34(s) + 192771 Vm4(s 



(B5) 



^W-o,)(^) ^ W|_J^(29p(,2(s) + 397nV(22(s) " 5pfi2(5) " HQpfU^) " 40p2l3(.^) + ^^^^^^^(s) + 5m'p^,,is) 



37r2 

+ 4^pfi2W) 

(B6) 



Pi'"^'\s) = y(99)'(p?io(^) - "^V(io(^)) (B7) 

pimiG.n^^^ = (99)(9Gg)(^pgo(,s) - ^mV(io(-) - ^™M2o(-^) - ^Pfi0(«) + %pU^) + §-Vfio(-) 

(B8) 

n(.-«)(«G.>2(^2 ) ^ 2^^-5)(5-G9)(^4n/(M2) _ n^/(M2)) (B9) 



10 



For the interpolating current 772: 
1 no r 



(BIO) 



Pi \^) — _4 



7 2 2 7 11 



57 1 

+ y™V(23(s) + 576TOVm(s) - 4to2p224(s) - 8toV'/33(s) - 384mVf44(s) - -^rn^ piizi^) 

+ 96m2(pf44(s) + pfi4(s)) - 48m2(p04(s) + pf44(s)) + 576mV(44(s) + ^pfi3(s) + ^"^Vf23(s) 

10 1 1 

- 96mVf34(s) - 96m4(pf44(s) + ^^14(5)) - ym^fisls) - 3P^;4(5) - 3^^224(5) 

/1 7 Q Q ^ 

- 120mVfi4(s) + ^Pn3(5) + 2"^Vf23(-') + yPm(s) + 80mVm('S) + ^^m^ p\^^{s) 

1 8 

+ 1536mVf45(s) + 2mV224(s) + ^Pi2a{^) + i^P21a{s) 



73 



^(,G,) (^^ = "^<gg^^^g) (iip(^^(,s) + pf^,(s) - 44p?2'3(.) + 8pi(3(.) + 9mV(22(s) " ypfi2(s) + 10mVf22(s) (g^2^ 
-"^Vfl2(s)) 



P^'"'^'(^)-|(9'Z)'(p?io(^)-"^Vio(^)) 



P2 



(qq)(qGq)l, 



17 



79 



17 



(B13) 



29 



29 



[s) = {qq){qGq){—p%^{s) - — mVno(s) - T^w^p{2o(s) " i;Piw{s) + j^Piwis) + -r^ni^ Piio{s) 



as 



36 



18 



n 



{qq)(qGq)2 



{Ml) = -{qq){qGq){m^Il\Ml) n"(A/|)) 



48' 



48 



(B14) 
(B15) 



For the interpolating current 773 
384 



pert ( \ 

P3 (s) 



TT 

a, 17 



(I6pfi5(s) + mVfi4(s) - 2m2pfi4(s) + 6mVu4(5) - P?i4(s) + 2pfi4(s) - 4p;^i4(s)) 
25 



+ -(-7-Pn5(s) - :tJ"iV(i4(s)) 



TT 6 



24 



(B16) 



P3 



(s'GG), X _ K9^GG) 



>{s) = 



23 



Vmls) + yPl23(s) - 3^213(5) - yP224(s) + gP?23(s) + 3P213(s) " yP224(s) 



10 



39 5 

- yP(23(s) + 1152m2p(34(s) - 20mVL4(s) " 16^^1^33(5) " 768m^p(4^{s) - -p'^^si^) 

+ 192m2(pf44(.) + pg,{s)) - 96m2(p04(.) + p^s)) + U52m*pi,,{s) - ^pfial^) + Q^'p^23i^) 

10 5 5 91 

- 192m4(pf44(s) + pfi4(s)) + ymVfi3(s) - 3P2''i4(s) - 3"^'pf24(s) - 240m2pf,4(s) + ^pfi3(s) 

- ^^[23(5) - yPm(s) + 160TOVm(s) + 3072mVM5(s) + Wm^P22iis) - y^P?24(s) - 3P?i4(s) 



4„L 



192m>f34(s) + 192mVi44(s 



(B17) 



piqGq) ^^^ ^rn{qga ■ Gq) 



3^2 (29p(i2(s) + 39mV(22(s) " ^^12(5) " l^GPiU^) " ^0p2i3is) + Wm'pi22{s) + ^^^P2i2{s) 



521 



Pn2(5)) 



(B18) 



11 



16 



pr (s) = y(99)'(prio(^) - "^v(io(«)) 



(B19) 



(B20) 



n 



(qq){qGq)2i,.r2 



(Ml) = -{qq){qGq){m^n\Ml) - n^^(A/|)) 



3 ---' 3 

For the interpolating current 774 

5 



pf (s) = ^[(I6pfi5(s) +™'pfi4(5) - 2m2pfi4(s) + 6m2p(i4(s) - p'^,^{s) + 2p^,M - Api,^{s)) 



^(^/5n5(5) - ^™Vn4(s)) 



(B21) 



(B22) 



pK^^>(.) = ^^ 



- -^nL. 



7 2 2 7 11 

3pf23(s) + gPl23(s) + 3^213(5) - 3^224(5) - gP?23(s) - gP^lsls) " 3^224(5) 

57 1 

-^™Vi23(s) + 576771^(34(5) - 47772^224(5) - ^rn^pissi^) - 334777^^44(5) - 3"^V223(5) 

+ 967n2(pf^^(s) + p^^^is)) - 487772(pO,(s) + pf,,{s)) + 576m'pi,,{s) + ^pf.^is) + ^mVf23(5) 

10 1 1 (^'^^ 

- 96777Vf34(5) - 967774(pf44(s) + p^,,{s)) - -,7l2pf^3(s) - -p^^is) - ^^^ P224{s) 

/1 7 Q Q ^ 

- 120777Vfl4(5) + ^Pn3(5) + 2"^Vm(5) + yPm(5) + 80777^^(5) + 96777 Vm4(5) 

1 Hi 

+ 15367«Vm5(5) + 2mV224(s) + -^P?24is) + 3/ofi4(s) 



^(.-G.>,.,_ 



*(5) = - 



m{qga ■ Gq) 



73 



37r2 
'^^^12(5)) 



(llp(i2(,s) + p^,^{s) ~ Upfi,is) + 8p',i,is) + 9m'pi,2{s) - -^P^^2is) + I0rn'pi22{s) 



(B24) 



pf'^ (5) = |(9?)^(p?io(5)-r77Vio(5)) 

(qq){gGq)l 
Pi 



(B25) 



(s) = (««)(«Gg)(— p/2o(5) - ^"^ Piio('S) - Ys"^ ^i2o(s) - 3^110(5) + ^/'iio('5) + ^"^ Piio(s)j 



n 



(qq) {qGq)2 , j.r2 



{Ml) ^ -{qq){qGq){m^U'{Ml) ~ n^^(A/|)) 



(B26) 
(B27) 



The functions p^^'-j; '"(s) and n^'^^(Af^) in the above expressions are defined as: 



pLkis) 



{-Ifi 



'"'"""' ^2r(/i)r(j)r(/fc)r(3 - /7 - j + fc) A^^„ '^"i/3_„ '^^ ^i+fc-7,^i+fe-j 



2-h~j+k 



-, (B28) 



^/^^^.■(■s) - 9 



(-1)'=4 



fe/l-fc-2 



7r2r(/7)r(j)r(fc)r(4 - /i - j + k) 



s2-h-j+k 



r/3™a. (a + /? - l)fe-i (7772(a + p) - apsY ^^ (2m'^(a + /3) + aps{h + j - k - 5)) 



(B29) 



12 



where h,j,k>0,h + j — k<2. 

I'hjkK'^} 7r2r(/i)r(j)r(fc)r(3 -h-j + k) 



l-h-j + k /„ 9, . n^ . n ,, . ■ , .^\ (B30) 



rP^.. {a + P- l)*^-! (to2(q, + ^3) - a^s)' " ^^^ (2m2(a + /3) + a/?s(/i + j-k- 4)) 



^i+fc-h^fc-j 



where /i, j, fc > 0, and h + j — k < 1. 



i (-l)*--4-'^-i p^ox /•'3-- (a + /3-l)fc-i(m2(a + /3)-a/3s)-''-J+'=^ 



[6(m2(a + ^) - al3s)^ - a/3s(6(77i2(a + /3) - a/3s)(2 + k - h - j) - a/3s{2 + k - h - j){l + k - h ~ j))] , 

(B31) 

where, h,j, fc > 0, and h + j — k < 0. 

^^^^^^) = 512.2r(4)2 I„„„ ^"X_ '^^ ^^ ^i_ ^" 36864(m2 - a.)^ ' 

(B32) 



""^"^ ' 7r2r(/i)r(j)r(A:)r(4 - /i - j + /c) 7„_ 7^_ "^ afe+i-Zi/^fc+i-j 

('2(6a - 1) (a/Js - m^{a + /3))^ - a/?s(3 - h - j + k) {a/3s{2a{h + j - fc - 8) + 1) + (12a - l)m^(a + /?))), 

(B33) 
where h,j,k>0, h + j — k < 1. 



^^'^^ ' 7r2r(/i)r(j)r(fc)r(4 - /i - j + fc) 7„„„„ 7/3„.„ a'^-'^/^fc - j + 1 

[6(8a - 3) {m^{a + /3) - a/3s)^ - 18(4a - l)a/3s(3 - h - j + k) {m^{a + P) - ajSsf (B34) 

+ 3(8a - l)a^li^s^{2 -h-j + fc)(3 - h - j + k) {m^{a + (3) - al3s) - 
2a^s^{l -h-j + fc)(2 -h-j + fc)(l - /i - j + fc)] , 
where h,j,k>Q, h + j — k < 1. 
^ , ^ to4 .a™„. a3s3(m2-(l-a)as)^ 



368647r2 7„^_ (m2 - asf 

1 /•"— /■'^— (a + /3-l)3((8a-3)TO4(a + /3)2--4a(10a-3);3m2s(a + /3) + 10a2(4a-l)/32s2) 

' da I dfi- 



20487r2r(4)2 J^^^^ J ^^^^ ' o?^ 

(B35) 

o 1 z-"™- a3s3(r„2_(^_^)Q,g^3 

P414(S) - 3ggg^ ^2 y^^^^ da ^2 (^2 „ ^5)3 

1 p™o. ^ /-/s^a. ^^ (a + /3 - \f ((8a - 3)m4(a + /3)2 - 4a(10a - 3)^m2s(a + /3) + 10a2(4a - l)/32s2) 

/3„ 



2048.2r(4)2 7„„^,„ ^^i. ^^ /34 



„ ,^ (_l)fc+i2-2fc-3 /-a,„ox /-/i^o. (a + /3-l)'=-i(m2(a + /3)-a/3s) '' ^'+'' ^ 



(B36) 



ijfew " ^2r(/i)r(j)r(fc)r(3-/i- j + fc)7„^^^^ 7^_ "^ a'=-''/3'=-j(i-a)2 

(^12(1 - a)^(8a - 1) {m^{a + 13)- afisf - 3 (48a^ - lOOa^ + 56a - 3) al3s{-h - j + k + 2) {afis - m^{a + I3)f 

-2{a- lf{24:a-l)a^l3^s^{h+j-k-2){h + j-k- 1) {aPs - m^{a + 13)) 

+ 4(a - l)^a'^P^s^{h + j - k - 2){h + j-k- l){h + j- k)) , 

(B37) 



13 



where h,j,k>Q, h + j — k < — 1. 



1 / 4m^ 

1 /S o\ / 4771,2 



^-(^) ^ -16^12 -'"^JV^-^' ^^^°) 



1 /S n\2 / 4771^ 



^m(^) = -Y^(2"'"') \/l-— ' (B41) 



1,1 / 4to2 



^-(^) - '8^( 4(1 -;/4.72) + v^ - "i^)' (^^') 



Pfio(^) - -^((3^ - 2-^)^1 - — - -TT^r^)' (B43) 

327r- ' V s ^s(s -4?772)' 

p , , -SSm*^ - 559m'*s + 326to2s2 _ 5353 

^"0 (*) ^ ion 2 r7 — r ^ ' (^^^) 

1207r2sY/s(s - 47n2) 

1 /"^ 771^ TJl^ 

n'(Ml) = j^X ■'";;;?M|°'pt- .(i-.)M| l- e^^^' 



The integration limits appears in the above expression are: 



_1 + ^1 -4m2/s _1- ^1 -4m2/s 

2 

/^mao; =1 - « /^mm = n ' (^^8) 

as — 771^ 



